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Abstract 

The magnitudes of the external gravitational perturbations associated with 
the normal modes of the Sun are evaluated to determine whether these solar 
oscillations could be observed with the proposed Laser Interferometer Space 
Antenna (LISA), a network of satellites designed to detect gravitational ra- 
diation. The modes of relevance to LISA — the I = 2, low-order p, f and 
g-modes — have not been conclusively observed to date. We find that the 
energy in these modes must be greater than about 10 30 ergs in order to be 
observable above the LISA detector noise. These mode energies are larger 
than generally expected, but are much smaller than the current observational 
upper limits. LISA may be confusion-limited at the relevant frequencies due 
to the galactic background from short-period white dwarf binaries. Present 
estimates of the number of these binaries would require the solar modes to 
have energies above about 10 33 ergs to be observable by LISA. 



1 



Typeset using REVTeX 



I. INTRODUCTION 



The oscillation modes of the Sun are now routinely observed by Doppler-shift mea- 
surements of the time-dependent velocity of the Sun's surface The frequencies and 
amplitudes of thousands of p-modes have now been accurately measured in this way. A 
number of detections of solar g-modes from Doppler-shift measurements have also been re- 
ported, but the claimed detections all have low signal-to-noise and they are not mutually 
consistent [[J. It has recently been claimed (although, again, it is not universally accepted) 
that certain periodic features observed in the solar wind provide a second, completely inde- 
pendent, means of observing the solar oscillations Some of the frequencies observed in 
the solar wind correspond closely to already-observed solar p-modes, while others are close 
to frequencies predicted for some solar g-modes. 

In this paper we investigate the possibility (also suggested independently by Schutz || 
and by Gough ||) that the solar oscillations might be observable using a third technique: 
measuring the external gravitational perturbations associated with these modes. The pro- 
posed Laser Interferometer Space Antenna (LISA) [7| is designed to detect the tiny gravi- 
tational fluctuations caused by distant sources of gravitational radiation with frequencies in 
the range 10 2 to 10 6 yuHz. This detector would also be capable of observing the fluctuations 
of the near-zone gravitational field produced by the normal modes of the Sun, if the ampli- 
tudes of the solar oscillations are sufficiently large. This new observational technique would 
be an interesting new probe of the structure of the Sun because it would be most sensitive 
for detecting a class of low-order modes that are presently unobserved using other meth- 
ods. In this paper we calculate the magnitudes of the external gravitational perturbations 
associated with the solar p, f, and g-modes to estimate whether these oscillations could be 
observed with LISA. 



II. SOLAR GRAVITATIONAL PERTURBATIONS 

We evaluate the coupling of the pulsation modes of the Sun to the external gravitational 
field using a reasonably accurate model for the structure of the Sun and its dynamics. 
We use the equation of state of a normal solar model (model 1 of Christensen-Dalsgaard, 
Proffitt, and Thompson ||) to determine the equilibrium structure of the Sun. We treat 
the pulsations as small amplitude perturbations which satisfy linear adiabatic evolution 
equations || with a realistic dynamical adiabatic index M. We solve the equations for 
these perturbations numerically to determine the magnitude of the external gravitational 
perturbation that results from a mode of given amplitude. The gravitational perturbation 
in the exterior of the Sun has the form 

' E \ 1/2 GM Q fRo\ l+ \, ^ 



\10 i8 ergs J R Q \ r J 

where r is the distance from the center of the Sun, and Y/ m is the standard spherical har- 
monic function. The factor a 2 8 representss the magnitude of this external gravitational 
perturbation for an oscillation mode normalized to have energy E = 10 28 ergs. The energy 
E is defined as 
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E = J p 5v*Jv a d 3 x, (2) 

where p is the mass density and 5v a is the (Eulerian) fluid velocity perturbation. The values 
of «28 that we compute for a number of solar oscillation modes are presented in Table |. 
The magnitude of the external gravitational perturbation of a mode scales as y/~E since E 
is quadratic in the perturbed fluid velocity. The values of a 2 s given in Table [I] illustrate 
that the largest surface gravitational perturbation occurs in the I = 2 (73-mode when modes 
excited to equal energy levels are compared. Also presented in Table | are the mode masses, 
M mo de, defined as the ratio of the mode energy E and the average surface velocity of the 
mode: 

E = J 5v a 5v* a sin 6d9d<f>. (3) 

The mode masses make it possible to convert the observable surface velocities of modes into 
mode energies. 

In our numerical solution of the pulsation equations we use a simple, but somewhat 
unrealistic, treatment of the surface of the sun. We set the (Lagrangian) perturbation in 
the pressure to zero at the point where the pressure vanishes in our solar model. For our 
equilibrium model we use the realistic equation of state || for densities above 10 _6 gm cm~ 3 . 
For densities below this value we smoothly attach a polytropic "atmosphere." The mode 
frequency uj and the gravitational parameter a^s are rather insensitive to the choice of 
density where this artificial atmosphere is attached. However the mode masses, M mo< j e , 
which depend on the surface values of the perturbed velocity, are rather more sensitive. We 
estimate that the errors in our computations of uj and a^s are about 0.01%, while the errors 
in M mode due to our simplified treatment of the surface may be as much as 10%. 



III. LISA 

LISA is designed to detect the presence of gravitational radiation, or other time depen- 
dent perturbations in the gravitational field, by monitoring the precise distance between 
pairs of satellites. Consider two satellites which form one arm of the detector. A near- 
zone gravitational perturbation (such as that produced by the solar oscillations) causes the 
distance between two satellites to oscillate with the frequency uj of the gravitational pertur- 
bation. The amplitude 5L of the periodic displacement of the satellites is determined by the 
standard tidal acceleration formula: 

uj 2 5L = Ln a n b V a V b 5®, (4) 

where L is the average separation of the satellites. The unit vector n a that appears in 
Eq. (D) gives the direction of the line between the satellites. To facilitate comparison with 
the published LISA sensitivity curves, we define the dimensionless strain h that would be 
measured by a detector consisting of two perpendicular interferrometer arms: 

h = uj- 2 {n a n b - m a m b )V a V b 5<5>. (5) 
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where n a and m a are unit vectors pointing in the directions of the two detector arms. The 
right side of Eq. @ has a complicated dependence on the orientation of the plane of the 
detector arms, the orientation of the detector arms within this plane, the position of the 
center-of-mass of the satellites, and the particular / and m of the mode. For simplicity we 
will use an average value for this quantity. We first average over all orientations of the two 
perpendicular arms, at a fixed point in space, to obtain 



< 



[n a n b - m a m b )V a V b 5<& 2 > = ^V a V 6 5$V a V 6 M>*. (6) 





We now average the right side of Eq. @ over all possible locations of the satellite around 
the Sun to find the mean-square value of h, 

< \h\ 2 > = — J V a V 6 5$V a V b 5$* smed9d<f). (7) 

Finally using Eq. (|T|) we evaluate the angular integrals on the right side of Eq. (0) to find 
the average strain h that would be sensed by the LISA detector 




2)(2/ + l)(2/ + 3) 
E \ 1/2 GM G fR Q ^ l+3 



10 28 ergs J uj 2 R% V r 



(8) 



The modes of the Sun with the same n and I but different m have frequencies that 
are slightly split by the Sun's rotation. In one year of integration, LISA will be able to 
distinguish frequencies to within about Af « 3 x lCT^ILz (for signal-to-noise of order 
unity). This Af is roughly an order of magnitude smaller than the size of the splitting due 
to the Sun's rotation. Thus modes with different m are effectively non-degenerate, and the 
energy E that appears in Eq. (0) must be that of a single mode of given I and m. 



IV. SIGNAL STRENGHS AND NOISE LEVELS 

To determine whether the oscillation modes of the Sun could be detected by LISA we 
must determine whether the expected gravitational signals exceed the expected noise in the 
detector and any "background" gravitational wave signals. We now give estimates for the 
parameters that determine the signal strengths and the noise levels relevant to this problem. 

As we have seen, Eq. (|5J), the magnitude of the signal in the LISA detector is determined 
by the energy contained in the modes of the Sun. The modes which are most likely to excite 
the LISA detector have not been observed to date. This makes the prospects of gravitational 
observations of these modes very interesting, but it also makes it very difficult to make 
reliable predictions of the energy confined in these modes. The observed low-/ p-modes have 
a maximum energy/mode of roughly 10 28 ergs for / near 3000 /iHz, with the energy/mode 
decreasing to roughly 10 27 ergs at / near 1000 /iHz ||10|| . A simple extrapolation of this curve 
would suggest even lower energies for the modes relevant to LISA (100 to 700 /iHz), but such 
an extrapolation might well be naive. Goldreich and Murray have argued that the decrease 
in p-mode energy (at fixed I) from 3000 to 1000 /zHz is due to the scattering of energy from 
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higher- n to lower-n p-modes of the same frequency [pjj . But the p-modes relevant to LISA 
are the lowest-n p-modes, so they would not be damped by this mechanism. 

At present only upper limits for the energies of these low-frequency modes are known. 
This is mostly because, for a given / and mode energy, the detection of modes by Doppler 
shift measurements becomes much more difficult as one goes to lower frequency modes. This 
difficulty is due to the fact that the the surface velocity (for fixed energy) decreases roughly 
like / 3 ' 5 for p-modes and f 2 for g-modes, while instrumental and background solar velocity 
noise both increase at lower frequencies ||. The present observational upper limit on the 
surface velocities of these modes is about 4 cm s _1 |3J. Using the mode masses computed 
in Table |, this translates into limits on the energies, -E max) which are also listed in Table | 
for each mode. It has been suggested |§ that the low signal-to-noise "detections" of solar 
g-modes are an indication that the g-mode energies are in fact just below these maximum 
values. These energies are rather large, and many Solar modes would be well above the 
LISA noise levels if their energies are close to these values. 

Kumar, Quataert, and Bahcall |12| have recently made predictions of the energies con- 



tained in the low-frequency solar g-modes. They consider the turbulent convection excita- 



tion mechanism used by Goldreich, Murray, and Kumar |L3| to explain the excitations of 
the observed p-modes. Applying this mechanism to the low-Z, low-n g-modes, they find that 
turbulent convection would excite surface velocities of about lCU 2 cm s^ 1 in these modes, 
assuming the dissipation time for the modes to be about 10 6 years. These surface velocites 
correspond (using our computed mode masses) to energies between 2 x 10 27 and 2 x 10 29 ergs 
for the g-modes listed in Table Q. Thus, a typical value for the energy predicted by this 
model is about 10 28 ergs. Similar energies are predicted for the low-/, low-n p-modes [|I~4| . 



For frequencies below 1000/iHz the anticipated noise in the LISA detector |15| would 
prevent the detection of periodic signals with amplitudes less than 



where / = u/2ir is the frequency of the signal ||16|| . It is also possible that other gravitational 
wave sources might produce a background "noise" from which the Solar oscillation signals 
could not easily be distinguished. Hils, Bender, and Webbink [[HJ have estimated that there 
are about 3 x 10 6 short-period white-dwarf binaries in our galaxy, producing a gravitational 
wave background of about 30 h sp in the relevant frequency band. At the time these estimates 
were made there were no known examples of white dwarf binaries with periods less that one 
day. Recently, however, two white dwarf binaries with periods of 3.47 and approximately 
4 hours resectively were discovered by Marsh and Marsh, Dhillon, and Duck [|19|| . Given 



these new observations, it should be possible to make considerably better estimates of the 
gravitational wave background from binaries, but to our knowledge, such estimates have not 
yet been completed. 



V. RESULTS 



Table | lists the values of the strain amplitude h, Eq. (||), for a number of modes of the 
Sun, assuming that each mode has an energy of 10 28 ergs and that the detector is located at 
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a distance r = 1 AU from the Sun. Only the / = 2 modes are considered, since the modes 
with successively higher I values have gravitational signals which are reduced below their 
I = 2 counterparts by additional factors of Rq/t »s 0.005. The modes in Table I are listed 
in order of increasing E min , the mode energy for which h = 3h sp . This is the minimum 
energy needed in a mode so that the mode would be observable above the detector noise 
with S/N = 3. Included in Table | is every I = 2 mode with E min < 10 32 ergs. We 
see that mode energies of at least 10 30 ergs are required for any of them to be observable 
by LISA. This is a factor of (up to) 10 3 smaller than the upper limits, but a factor of 
10 2 larger than the energies estimated on theoretical grounds by Kumar, Quataert, and 
Bahcall [jE^] . Given these energies, the corresponding surface velocities of these modes would 
be about 10 _1 cm s -1 , which is roughly the detection threshold for the recently-launched 



SOHO satellite pD |. This indiciates that any modes that could be observed by LISA would 
already have been observed by SOHO. Nevertheless, LISA could still provide some unique 
information since it measures the gravitational amplitude of the mode which depends on the 
density perturbations throughout the interior of the Sun. 

If the background radiation due to short-period white dwarf binaries is comparable to 
the values predicted by Hils, Bender, and Webbink Ui~7| , then at / ~ 300 /zHz the gravita- 
tional wave background is a factor of about 30 larger than LISA's detector noise. Hence 
mode energies about 10 3 times greater than the -E m i n 's listed in Table | would be required for 
solar oscillations to rise above this background. Such energies would be close to the present 
upper limits for these modes. Therefore, if the gravitational wave background is close to the 



strength predicted by Hils, Bender and Webbink [[17]], then the mode energies would have 



to be very close to their maximum allowed values for LISA to contribute at all to helioseis- 
mology. LISA's one-year orbit around the Sun would modulate the solar oscillation signal 
in a way that differs from the one-year modulation of the background gravitational wave 
signal. It is conceivable that this difference could be used to help distinguish solar oscilla- 
tions from the gravitational wave background, but we have not investigated this possibility 
in any detail. 
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TABLES 



TABLE I. Gravitational Parameters of Solar Oscillations. 



Mode 


f (A*Hz) 




«28 




M modc (gm) 




h 




^min( 


erg) 


^max 


(erg) 


93 


220.4 


7.17 


X 


io- 


-12 


8.3 x 10 31 


1.18 


X 


10" 


-23 


1.55 x 


10 30 


1.3 x 


10 33 


92 


254.0 


6.07 


X 


10" 


-12 


4.4 x 10 31 


7.50 


X 


10' 


-24 


1.97 x 


10 30 


7.0 x 


10 32 


Pi 


381.6 


5.31 


X 


10" 


-12 


4.3 x 10 30 


2.91 


X 


10' 


-24 


2.02 x 


10 30 


6.9 x 


10 31 


9a 


192.2 


5.96 


X 


10" 


-12 


1.9 x 10 32 


1.29 


X 


10" 


-23 


2.42 x 


10 30 


3.0 x 


10 33 


95 


168.8 


4.72 


X 


10" 


-12 


3.5 x 10 32 


1.32 


X 


10" 


-23 


4.19 x 


10 30 


5.6 x 


10 33 


9i 


293.6 


2.81 


X 


10" 


-12 


2.3 x 10 31 


2.60 


X 


10" 


-24 


8.44 x 


10 30 


3.7 x 


10 32 


f 


350.9 


2.31 


X 


10" 


-12 


1.2 x 10 31 


1.50 


X 


10" 


-24 


1.12 x 


10 31 


1.9 x 


10 32 


P2 


514.4 


2.02 


X 


10" 


-12 


5.5 x 10 29 


6.07 


X 


10' 


-25 


1.17 x 


10 31 


8.8 x 


10 30 


96 


134.0 


2.80 


X 


10" 


-12 


8.2 x 10 32 


1.24 


X 


10" 


-23 


1.37 x 


10 31 


1.3 x 


10 34 


P3 


663.6 


8.60 


X 


10" 


-13 


1.4 x 10 29 


1.56 


X 


10" 


-25 


5.53 x 


10 31 


2.2 x 


10 30 


97 


101.2 


1.42 


X 


10" 


-12 


1.6 x 10 33 


1.10 


X 


10" 


-23 


6.30 x 


10 31 


2.6 x 


10 34 
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